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W) 

The N cardinality k ideals of any w-element poset (w, k variable) can be enumerated 
in time 0(Nw 4 ). The corresponding bound for fc-element subtrees of a w-element tree is 

y—{ 0{NW 5 ). 

o 

O 1 Introduction 

■i— > 

a 

An implication is a Boolean formula of type 

{ (1) (ai A a 2 A • • • A a m ) — > (b\ A b 2 A • • • A b n ). 

> 

Each finite closure system C can be viewed as the set of models of a suitable family E of 
implications, called an implicational base. It was known [KS] that from £ one can enumerate C 
£N1 in output-polynomial time; see [W2] for a "concrete" algorithm. Here we show that for any fixed 

OO k the A:-element members of C can be enumerated in output-polynomial time as well, provided 

C satisfies some technical condition (Theorem 2). As Corollary 3 and Corollary 4 we obtain the 
two results mentioned in the abstract. 



• *-h Before we come to the section break up, let us give some background on Corollary 3. For an 

arbitrary ^-element poset (W, <) we shall write N for the number of its (order) ideals, and Nf- for 
the number of /c-element ideals (0 < k < w). Calculating iV is a #P-complete problem, but for 
special cases, such as 2-dimensional posets or interval orders, one can do it in u>-polynomial time. 
As to enumerating (= producing) all N ideals, obviously this can't be done in ^-polynomial 
time since for general posets (e.g. antichains) N grows exponentially with w. However, by the 
mentioned [KS] result it can be done in output-polynomia¥\ time. In fact Medina and Nourine 
[MN] enumerate the ideals in time O(Nd) where d is the maximum number of covers of an 
element of W . 



Let us turn to ideals of fixed cardinality k. Calculating Nk can be done in it;-polynomial time 
for interval orders [S], and for some posets explicite closed formulae for iV^ are given in [C]. As 
to enumerating the ideals of cardinality k, for interval orders it was known to be doable 

"Thus in time 0{N c w d ) for some c, d £ Z+, as opposed to say 0(A rv/ ™) which is not output-polynomial. An 
output-polynomial algorithm which has c = 1 is called output-linear. 
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in ./Vfc-linear time [HNS]. We show in Corollary 3 that output- linear time can be achieved for 
general posets. 

Section 2 introduces, on a toy example, the 0, 1, 2-algorithm which underlies the proof of Theorem 
2 in Section 3. The 0, 1, 2-algorithm is convenient for that particular proof, but gets fine-tuned 
in Section 4 to the (a, 6)-algorithm which targets only singleton premise implications (so m = 1). 
This kind of implication, which is intimately linked to posets, is the sole kind we are concerned 
with in the remainder of the article. The number of A:-element ideals of a poset (W, <) can 
be viewed as the fc-th Whitney number of the associated distributive lattice D. In Section 5 we 
pit the (a, 6)-algorithm against a straightforward recursive method to get the Whitney numbers. 
Section 6 indicates variations of the (a, 6)-algorithm, for instance whether it carries over from 
implications like a — > {b\ A&2A&3) to other kinds of "quasi-implications" such as a — > {b\ A62A&3). 



2 The 0, 1, 2-algorithm in a special case 

Rather than (1) we shall mostly adopt the set theoretic notation {ai, • • • , a m } — > {b% , ■ • • , b n } 
for implications. Let us introduce the 0, 1, 2-implication algorithm with an example where all 
implications have m = 1. The adaption to m > 1 happens in Section 3. 

In order to calculate the family Mod of all /c-element ideals of a poset (W, <) , we represent subsets 
of W by their characteristic bitstrings of length w but with the additional proviso of a don't 
care symbol 2. Thus the powerset V(W) is (2, 2, • • • , 2), and a general {0, 1, 2}-valued row, like 
r = (1,2,2,0,1,2,0,2,1), can be viewed as an interval in V{W): Taking W = [9] := {1, 2, • • • , 9} 
we have r = {X G V{W) : {1,5,9} C X C {1,2,3,5,6,8,9}}. Consider the poset (W, <) in 
Fig.l. with W = [6]: 



Obviously the order ideals of (W, <) are exactly the subsets X C W that simultaneously satisfy 
all implications 



Figure 1 
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(2) 



{3}^{1}, {4}->{l,2}, {6} ^{3,4}, {5} ^{4} 
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in the obvious sense that if an implication {a} — > B happens to have its premise {a} in X, then 
X must also contain its conclusion B. The upcoming 0, 1,2- algorithm adheres to the principle 
of exclusion [W2] in that starting with the powerset (2, 2, •••,2) the implications (generally 
: constraints) are "imposed" one by one until Mod arises as a disjoint union of 0,1, 2- valued 
(generally: multivalued) rows. 
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Table 1 



After imposing the first implication {3} — > {1} on (2,2,2,2,2,2) the working stack consists of 
n = (2,2,0,2,2,2) and r 2 = (1,2,1,2,2,2) because n contains all X's with {3} £ X, and r 2 
contains all X's with {1,3} C X. Note that r\ and r 2 encode disjoint set systems. For both 
rows the next pending implication is the second one (PI = 2), but always only the top row 
of the working stack is handled: This is known as last in, first out (LIFO) stack management. 
Imposing {4} — > {1,2} upon n yields the sons (2,2,0,0,2,2) and (1,1,0,1,2,2), and so forth. 
As soon as a row is finalized, i.e. the last implication has been imposed on it, it is put on the final 
stack. Here this happens the first time when {5} — > {3} is imposed on (2, 2, 0, 0, 2, 0), making it 
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the finalized row (2, 2, 0, 0, 0, 0). We leave it to the reader imposing the third implication to the 
new top row (1, 1, 0, 1, 2, 2), and to carry on a few steps until one has the working stack empty 
and the final stack comprising four rows. Due to the workings of the algorithm, the rows are 
disjoint, and so we can add up their cardinalities and get 

(3) |Mod| = 4 + 2 + 2 + 4 = 12 

For later use we state the following. 

Lemma 1: Let r\ be a {0, 1, 2}-valued row of length w, and let k E [w] be fixed. Then the M\ 
members X £ r\ with |X| = k can be listed, given that M\ > 1, in ordinary set notation in time 
0{M lW ). 

Proof. Say w = 8 and n = (0,0,2,2,1,2,0,2) and k = 3. Putting s = |ones(n)| and t = 
|twos(ri)| yields M\ = ^ * ^ = , namely these six sets: 

X = {5}UY (Y C {3, 4, 6, 8}, \Y\ = 2). 

From this it is clear that for general r\ the M\ sets X G r±, with \X\ = k can be listed in time 
0{M lW ). ■ 



3 The flag of Papua and the main result 



The algorithm of Section 2 can be enhanced as follows to handle arbitrary implications (m > 1). 
Say A — > B is {1, 2, 3, 4, 5, 6} — > {8, 9} and it is to be imposed upon row r in Table 2. If one of 
its symbols * at the first or second position is then all X G r satisfy A — > B. This is also the 
case if both the eighth and ninth * is 1. Let us assume none of that occurs. Then not all sets 
in r satisfy A — > B and the ones that do belong to either 

S~ :={Xer : A%X} or S + := {X G r : A U B C X}. 

A moment's thought confirms that 5_ is the disjoint union of candidate sons r\ to r±. 
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Table 2 



As to £+, if either the 8 th or 9 th component of r is then S + = 0. Otherwise S+ is modelled 
by the loner candidate son rs (those among the 1 st , 2 nd , 8 th and 9 th components of r that are 1 
stay that way in r$, and those which are 2 turn to 1). 
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Generally, if an implication A — > B is to be imposed upon a row r which has zeros(r) n A = 
and |twos(r) D A\ = t then the boldface pattern in Table 2 becomes a t x t matrix that has 
l's below and 2's above its 0-diagonal. That's akin to the flag of Papua (New Guinea), which 
provides a shorthand symbol for the situation. 

For a family T of closure systems C with universes of unbounded cardinality w consider the 
following disjoint extension property: 

(4) There is c > 1 such that for all C G J 7 , all Z G C, all Y C W and all k < w, it can be 
checked in time 0(w c ) whether there is Z G C with Zq C Z and \Z\ = k and F PI Z = 0. 

(If already Y n Zo / or if \Zq\ > k, the answer is easy.) 

The length of an implication A — > B is \A\ + and for a family S of implications ||S|| is the 
sum of the lengths of its members. Observe that ||S|| < Further let s(S) < w be the 

maximum cardinality \A\ occuring. 

By [RM, Thm.10.3] computing the closure Zq of the empty set costs 0(||£|| + w). By (4) 
(putting Y = 0) it costs another 0(w c ) to decide whether C contains any fc-element member Z. 



Theorem 2: Let J 7 be a family of closure systems that enjoys property (4). Then for each 
C G T that comes with an implicational base S and each k G [w] the M sets IgC with 
\X\ = k can be enumerated in time 0(Mti;s(S)|S| 2 + Mw c+1 \Y,\). 



Proof. Let C G T be based on W = [w] and fix k G [w;]. Also fix an arbitrary ordering 
Ai ->■ Bi, A 2 -)• -B 2 , ■ ■ ■ of the implications in S. Starting with r = (2, 2, • • • , 2) = V{W) we 
always focus on the top row r of our working stack and impose its pending implication (PI), 
very much akin to the previous example. 

Specifically, if there isn't X G C with \X\ = k, this can be discovered in time 0(||S|| + w c ) as 
seen and we output M = 0. Otherwise by induction rooted in row (2, • • • , 2) we can assume 
that the working stack's top row r is feasible in the sense of containing at least one Z G C 
with \Z\ = k. Using a threefold case distinction we next show that upon imposing the pending 
implication A — > B (we dropped its index) on r, at least one of the arising candidate sons will 
remain feasible. 

Case 1: Either A n zeros (r) ^ or B C ones(r). Then all X G r satisfy A — > B. Thus r 
remains unaltered but has its Pi-pointer raised by 1. 

Case 2: A n zeros(r) = and B n zeros(r) / 0. Then the sets X G r satisfying A — >■ P» are 
exactly the ones in S* _ = {X G r : i ^ X}. In particular Z G 6*" (since Z G C entails that Z 
satisfies ^4 — > B). Hence if r\, ■ ■ ■ , rt (t < s = s(S)) are the Papua candidate sons of r triggered 
by S~ then at least one r\ is feasible (i.e. the one with Z G r« and perhaps others). As to 
checking feasibility, see (5) below. 

Case 3: A n zeros(r) = and B C ones(r) U twos(r) but 1? £ ones(r). Then the sets X G r 
satisfying A ^ B are exactly the ones in = {A G r : A % X] or S" + = {A G r : AUB C I}. 
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Here both sets are nonempty. Apart from the t Papua candidate sons there will be a loner 
candidate son encoding S + with l's written all over AU B. Again at least one of these t + 1 
rows is feasible. 

Here comes the core claim: 

(5) For any 0, 1, 2-valued row f it can be decided in time 0(||S|| + w c ) whether or not it is 
feasible. 

Proof of (5). Again by [RM] the C-closure Zq of ones(r) C W can be computed in time 0{\ |E| | + 
w). Each closed set contained in f comprises ones(r) and hence comprises Zq. Thus if \Zo\ > k 
then r is infeasible. If \Zq\ < k put Y = zeros(r). By assumption one can test in time 0(w c ) 
whether or not there is a closed /c-element set Z that comprises Zq and is disjoint from Y (and 
thus Z G r). Only in this case is r feasible. All of that costs 0(||S|| + w c ). 

We come to the overall cost assessment. Let R be the number of final rows. Because they are 
disjoint, R < M. Since deletion of intermediate rows is avoided, the total number of impositions 
occuring is at most i?|£|. The cost of one imposition is as follows. Making the t or t + 1 
candidate sons r-i in Case 2 or 3 costs 0((t + l)w) = O(sw), and checking each for feasibility 
costs 0(s||S|| + w c+1 ) by (5). 

Each final row n solely consists of closed sets. By Lemma 1 sieving the M, > 1 members X G r\ 
with \X\ = k costs O(Miw). Doing this for all R final rows costs 0(M\w) + • • • + O(Mrw) = 
O(Mw). Adding to it the cost of all impositions yields 

0{Mw) + 0{R\T\ ■ (s||S|| +w c+l )) = 0(Msw\Z\ 2 + Mw c+1 \Z\). 

□ 

Corollary 3: For any fixed k G [w] the M = many /c-element ideals of any w-element 
poset (W, <) can be enumerated in time 0{Mw 4 ). 



Proof. An implicational base S for the closure system C of all ideals of (W, <) consists of all 
implications p — > LC(p) where p £ W is non- minimal, i.e. has a nonempty set LC(p) of lower 
covers. Finding S costs 0(w 2 ). 

Let Zq be an ideal of (W, <) and let Y C W be disjoint from Zq. If Y' is the filter generated by 
Y (which remains disjoint to Zq) then each ideal disjoint from Y is contained in W \ Y' . Put 
P = W\(ZqL)Y'). If \Zq \ + \P\ < k there is no Z as required. Otherwise we can shell from below 
any k — \Zq\ elements of P and join them to Zq. This yields an ideal Z as required. Calculating 
\P\ costs 0(w 2 ), and so the claim follows from Theorem 2 in view of s(S) = 2, c = 2, |S| < w. 
□ 

It is an open question whether "fc-element ideals" in Corollary 3 can be generalized to "ideals 
of weight fc" with respect to some weight function W — > Z + . If the answer is affirmative, 
that would entail (by considering the factor poset of strong components) that for any family 
S of singleton-premise implications (possibly with directed cycles) the S-closed /c-element sets 
can be enumerated in output-polynomial time. We mention that finding one maximum weight 
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/c-element ideal of a poset is NP-hard [FK] . 



Corollary 4: Given a ^-vertex tree T and k G [u;], the M = many /c-subtrees of T can be listed 
in time 0(Mw 5 ). 



Proof. An implicational basis S for the closure system C of all (vertex sets of) subtrees of T 
consists of all implications {01,02} — > {bi,--- ,b n } where a±,a2 are non-adjacent vertices of T 
and a±, b\, ■ ■ ■ , bt, a<i (t > 1) is the unique path from a\ to 02- One can find S in time 0(w 2 ). 

Let now Zo be a subtree of T and let F C V(T) be disjoint from V(Zq). Let F be the subforest 
of T with = V(T) \ Y. There is a subtree Z of T with V(Z ) C V{Z) and |F(Z)| = k if 

and only if the connected component of F containing X has cardinality > k. Finding the size 
of the connected component of a vertex in a graph can be done in linear time with respect to 
the number of edges. In our case that's 0(w), and so the claim follows from Theorem 2 in view 
of a(E) = 3,c = 1,|E| < w 2 . □ 

There is some literature on counting or enumerating various types of binary trees and their 
subtrees. As for arbitrary trees T, all subtrees of T can be enumerated in output-polynomial 
time [R] but Corollary 4 seems to be new. Observe that different from posets a maximum weight 
A;-subtree of T can be found [FK] in time 0(u> 4 ). 



4 The (a, 6)-algorithm 



Suppose the implication {1, 2, 3} — > {5} is to be imposed upon (2, 2, 2, 2, 2). Instead of the three 
Papua candidate sons and the loner, we could just take (n, n, n, 2, 2) and the loner (1, 1, 1, 2, 1). 
Here the wildcard (n, n, n) by definition means "at least one here" . The resulting implication 
n-algorithm [W2] indeed beats the 0, 1, 2-algorithm in practise but for the proof of Theorem 2 the 
0, 1, 2-algorithm was more convenient. For singleton-premise implications the trick (n, ■ ■ • , n) 
(1, • • • , 1) boils down to e 1, but there is another way to improve upon the 0, 1, 2-algorithm in 
this case. In brief, we shall employ two new symbols a, b. Rather than imposing say {3} — > {1, 5} 
on (2, 2, 2, 2, 2) by splitting it in two, we replace (2, 2, 2, 2, 2) by a single row (6, 2, a, 2, b) which 
by definition represents the intended set system. 

Before we embark on the ensuing algorithm, we formally define an {0, 1, 2, a, b} -valued row r as 
a partition of [w] into four (possibly empty) parts zeros(r), ones(r), twos(r), implications(r), 
such that if implications(r) 7^ (f>, it is further partitioned as 

implications (r) = premfl] U conc[l] U • • • U premfs] U conc[s] (s > 1), 

where all prem[i] are singletons, and all conc[i] are nonvoid (1 < i < s). A set X C [w] by 
defintion belongs to r if and only if ones(r) C X and zeros(r) n X = 0, and for all 1 < i < s it 
holds that 

prem[i] € X =>• concfi] C X. 

While this defintion reflects the author's Mathematica-implementation of {0, 1, 2, a, 6}-valued 
rows in the upcoming a, 6-algorithm, we shall use a more visual representation. Thus, up to 
permutation of the entries, a typical {0, 1, 2, a, 6}-valued row looks as follows: 
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(6) r = (0,0, 1,1,2, 2,2, ai, &i, h, a 2 ,b 2 ,b 2 , b 2 , a 3 ,6 3 ). 

For instance, ones(r) = {3,4} and prem(2) = {11}, conc(2) = {12,13,14}. It is clear that 
generally 

(7) \r\ = 2l twos ( r )l • f[(l + 2l prcm «l) [where J]( ■■) = !)■ 

i=i V i=i / 

Before we tackle the details of the a, 6-algorithm in the proof of Theorem 5 let's glimpse how it 
processes the implications in (2). Table 3 is quite self-explanatory: 
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Table 3 

Notice that the numbers of rows both in the working and final stack have shrunk considerably 
compared to Section 2. Using (7) one confirms (3), i.e. that 

|Mod| = 2 -(1 + 2) + 2 -(1 + 2) = 12. 



Theorem 5: The (a, 6)-algorithm generates all N order ideals of a w-element poset (W, <) as 
a disjoint union of R < N many {0, 1, 2, a, 6}-valued rows in time 0(Rw 2 ). 



Proof. We may assume that W = [w] and that 1, 2, • • • , w is a linear extension of (W, <) (we 
didn't have that comfort in the proof of Theorem 2). Fix a £ W with corresponding pending 
implication {a} — > B that needs to be imposed on some row r. In view of the linear extension 
ordering, r must have a 2 at position a and B is contained in {1, 2, • • • , a — 1}. 
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First case: There is a within the range of B. Then the 2 at position a becomes 0. 

Second case: All the S-positions are 1. Then r carries over unaltered and, as for all successful 
impositions, PI = a is raised to PI = a + 1. 

Third case: At least one 5-position is 2, and the others, if any, are 1. Then the 2 at position a 
becomes the symbol aj, and the 2's within B become symbols bj (where j is a fresh index). 

Fourth case: There are no 0's but symbols Oj or bi within the range of B (distinct implications 
having distinct subscripts). A typical situation would be: 

r =(... &! 6 2 , 0, 64,62,01,02,63,2, 63, 03, 04, 2 2, 2) 

V v ' 

B a 

r' = (•■■, 61, 6 2 , 0, 64,62,01,02,63,2, 63, a 3 , a 4 , 0, 2, 2) 
r" = (-.., 1, 1, 0, 1, 1, 1, 1, 1, 1, 63, a 3 , 2, 1, 2, 2) 

One checks that the sets X G r that satisfy {a} — )• S are precisely the ones in the (disjoint) 
union of r' and r" . As to r", notice that putting one 63 = 1 in 036363 results in 0363, whereas 
putting 64 = 1 in 0464 results in 2. 

Since the deletion of rows is avoided, the total number of impositions occuring is at most Rw 
as in the proof of Theorem 2. The cost of all impositions hence is 0(Rw 2 ) since imposing one 
implication (even in the fourth case) costs 0(w). The height of the LIFO working stack is easily 
seen (as in [W2, Thm.6]) to stay below 2w + 1 throughout. □ 

Enumerating the ideals of a poset has many applications in scheduling theory and other parts of 
operations research [CLM]. As to applications in pure algebra, one concerns the calculation of 
all submodules of a finite .R-module from its join-irreducible submodules. In fact that extends to 
calculating arbitrary finite modular lattices. Here one benefits from the flexibility of the (a, 6)- 
algorithm to target only specific "closed" ideals. Details will be elaborated elsewhere. Another 
application [W3] is the production of all pairwise nonisomorphic incidence algebras that share 
some common scaffolding. 



5 Calculating Whitney numbers of distributive lattices 

The number Nk of rank k elements of any ranked lattice L is called the fc-th Whitney number 
of L. Here we focus on the distributive lattice D(P) of all order ideals X of a ^-element poset 
W = (W, <), in which case the rank of X equals its cardinality. Thus Nk can be calculated in 
time 0(NkW 4 ) (Corollary 3). We refer to the introduction for various special cases concerning 
Nk- In practise the (a, 6)-algorithm beats the 0, 1, 2-algorithm underlying Corollary 3. Even 
more so when the issue is not enumeration but merely the calculation of N^. But in order to 
do so the (o, 6)-algorithm needs a little update. Namely, given a {0, 1, 2, o, 6}-valued row r, we 
wish efficiently to calculate 

(8) Card(r,A;) := \{X G r : \X\ = k}\ 
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for all k £ [w]. Associate with r a polynomial pol(r, x) with indeterminate x that comes as a 
product as follows. Each symbol 1 contributes a factor x, and each symbol 2 a factor 1 + x. 
Further, each pattern aibibi • • - hi with t = |prem(z')| contributes a factor 

i+(;>+(<y+- + (;>'+*■« 

It is easy to verify that Card(r, k) is the coefficient of x k in pol(r, x). For instance, r in (6) has 
pol(r,x) = x 2 (l + x) 3 (l + 2x + x 2 + x 3 )(l + 3x + 3x 2 + x 3 + x 4 )(l + x + x 2 ) 

= x 2 + 9x 3 + 37x 4 + 93x 5 H + 6x 13 + x 14 

and thus e.g. Card(r, 5) = 93. 

The multiplying out of factors is achieved "fast" using the Mathematica command Expand 
[■ • •] and we omit to delve into its theoretical complexity here (for some e-algorithm handling 
hypergraph transversals this was done elsewhere). Summing up the polynomials of all final rows 
yields the rank polynomiaf] 

w 

R(x,D) = J2 Nk °° k 

k=0 

of our distributive lattice D = D(W). 

Let us present another, this time recursive way to calculate R(x) = R(x,D). For any a G W 
put T_ = {X 6 D : a X} and T + = {X G D : a € X}. Because T_ equals D(W \ a f), and 
X i — y X \ (cb is a bijection from T+ onto D(P \ a \.), it follows^] that 

(9) R(x) = R-(x) + x a R + (x) 

where R-(x) and R + (x) are the rank polynomials of D(W \ a t) and D(W \ a |) respectively, 
and a = \ a I \ . 

Iterating (9) on the smaller posets W\a 4- and W\a t eventually yields R(x) explicitely. Here it 
pays to always pick a in such a way that \a ^ \ + \a t | is maximum. Furthermore, for n-element 
antichains we can dispense with (9) since their rank polynomials are obvious. (That's more 
generally true for disjoint unions of chains but they are hard to identify.) 

In Table 4 below we pit the (a, 6)-algorithm against the recursive method on some random 
examples. Specifically the lattices D are the ideal lattices of posets W(m, £, s) which consist of 
t levels Lev(i), all of cardinality m, in such a way that each a £ Lev(i) (2 < i < I) is assigned 
s random lower covers in Lev(i — 1). ■ 

Both algorithms were implemented with Mathematica 8.0, the implementation of the recursive 
method being straightforwarcPl In all cases both methods agreed upon R(x, D), thus very likely 



^This is not to be confused with the order polynomial of (W, <) whose coefficients count the number of 
homomorphisms from W to chains. 

*The recursion (9) has likely been discovered before but the author couldn't pinpoint a reference. We mention 
that a similar recurrence used in [BK, p. 106] to count the antichains of (W, <), in fact is the standard recurrence 
to calculate the independence polynomial of a graph (in this case the comparability graph of (W,<). Notice 
that while antichains and order ideals of (W, <) are in bijection, this is not the case for the respective fc-element 
objects. 

^ It is harder to set up and implement recursions for counting the fc-element subtrees of a tree; cf. Corollary 4. 
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both are correct. Rather than R(x,D) we record the total number N = R(1,D) of ideals, the 
respective times, as well as the more telling parameters nfr and nsum. The former is the number 
of final rows generated by the (a, 6)-algorithm, the latter the number of final, (i.e. coupled to 
antichains) summands constituting R(x, D) when using the recursive method. Depending on 
the type of poset one or the other algorithm prevails. Recall that the (a, 6)-algorithm has the 
additional benefit of delivering the ideals themselves in a compact format. 



(m, £, s) 


N 


time- (a, b) 


nfr 


nsum 


time-rec 


(5,28,2) 


30416508 


25 


8932 


4827 


7 


(3,28,1) 


66506610198 


61 


56256 


595830 


1070 


(10,15,9) 


15560 


3 


984 


228 


1 


(15,10,6) 


709150 


80 


25302 


5870 


6 


(20,2,2) 


500466940 


3 


1872 


30205 


49 


(20,3,2) 


27003561664 


746 


1041866 


946439 


1612 


(20,4,2) 


214942410560 


1335 


1148952 


6713622 


11892 



Table 4 



6 About quasi-implications 

Each singleton-premise implication, say {a} —> {b±, 62, 63} properly speaking is the propositional 
formula in (i) below. It serves as point of departure in our attempt to carry over the a, b- 
algorithm to other kinds of (singleton-premise) quasi-implications. We investigate the eight 
variations obtained from (i) by independently negating the premise a, respectively negating all 
variables 6i,i>2>&3 of the conclusion, respectively switching A to V. Since negating all variables 
can't deliver anything logically new (yet psycho- logically) , we mainly focus on (i) to (iv) rather 
than (i*) to (iv*): 



(i) 


a — 


> (61 A b\ A 63) 


(1,1) implication 


(ii) 


a — 


► (61 V6 2 V6 3 ) 


(0, 0') quasi-implication 


(iii) 


a — 


► (61 V b 2 V 62) 


(1,0') quasi-implication 


(iv) 


a — 


>■ (61 A T>2 A 63) 


(1,0) quasi-implication 


(i*) 


a — 


> (61 A b~2 A 63) 


(0, 0) quasi-implication 


(ii*) 


a — 


>• (61 V 6 2 V 63) 


(1,1') quasi-implication 


(iii*) 


a — 


>• (61 V 6 2 V 63) 


(0, 1') quasi-implication 


(iv*) 


a — 


> (bi A b 2 A 63) 


(0, 1) quasi-implication 



As to (ii), this is equivalent to (b\ A 62 A 63) — > a. Such singleton- conclusion implications are more 
powerful than singleton-premise implications since each implication, say (61 A 62 A 63) — > (aiAa2), 
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is equivalent to a conjunction of such implications (fei Ab 2 A 63) — > a\ and (61 Ab 2 A 63) — > 02- 
The (a, 6)-symbolism does nolj^] carry over to implications with premises of cardinality > 1. 

As to (hi), this is equivalent to the clause a V 61 V b 2 V 63. Also (ii) could have been written as 
a V61 V62 V63. By definition a conjunction of clauses with at most one positive literal, is called 
a Horn formula. The Horn n-algorithm of [W2] extends the implication n-algorithm in obvious 
ways. Having at most one negative literal per clause yields Anti-Horn formulae (cases (ii*) and 
(hi*)). 

As to (iv), consider the task to determine all independent sets of vertices (anticliques) in a 
graph G. Obviously X being an anticlique amounts to the satisfisfaction of all clauses a V b 
where {a, b} ranges over the edges of G. Grouping together all edges incident with a fixed vertex 
a this becomes (say) 

(a V 61) A (a V b 2 ) A (a V 63) = a V (61 A b 2 A 63) = a ->■ (61 A b 2 A 63) 

In [W4] this gives rise to the (a, b)- algorithm which bears a superficial resemblence to the (a, 6)- 
algorithm. However, (i) and (iv) are not duals in any sense, and also the technical details 
differ considerably. For instance, there is no such thing as a linear extension for graphs. As to 
(iv*), albeit logically and algorithmically equivalent to (iv), it e.g. handles this scenario: When 
composing a committee, the lack of person a (some VIP) necessitates that at least persons 
bi, b 2 , 63 must make up to that. 

Apart from (i) to (iv) versus (i*) to (iv*), one can distinguish the quasi-implications that feature 
A in their conclusions of fixed type (carrying tags (x, y)) from the ones that feature V (carrying 
tags (x,y')). Acyclic families of quasi-implications invite conglomeration to rooted trees, and 
this has indeed algorithmic benefits for type (x, y). Let us sketch the basic idea for type (1, 1). 
The conjunction 

(10) (1-)- (2A3A4)) A (2 — >■ (5 A 6)) A (4 -»• 7) A (5 (8 A 9)) 

of four implications gives rise to the (rooted) tree in Figure 2(a). 

^The attempt to carry it over degenerated in a messy case distinction and thus spawned the implication 
n-algorithm touched upon at the beginning of Section 4. 
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Figure 2 (a) 



Figure 2 (b) 
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We wish to encode with appropriate symbols * the set of all X G 7- > ([9]) that satisfy (10) 
as just one "tree-valued" row r = (*,*,*,*,*,*,*,*,*). Notice that {X G r : 1 G X} = 
(1,1,1,1,1,1,1,1,1) whereas {X G r : 1 X} = (0, *, 2, a, *, *, b, *, *) (see Figure 2(b)). 
Details pending, the aim is to impose trees T upon tree-valued rows r' by splitting r' in 
and r' 2 according to whether T's root is turned to 1 or 0. The resulting "tree-algorithm" would 
comprise our (a, 6)-algorithm as the special case where all occuring trees have height 1. Another 
direction worthwile exploring is the exploitation of potential symmetries present in E. 

We indicated algorithms for handling all eight types individually, but what if some types occur 
simultaneously? A case in point is [Wl] where (i), (i*), (iv), (iv*) do just that. 
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